Synthetic Spin-Orbit Coupling in Two-level Cold Atoms 
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Synthetic spin-orbit coupling (SOC) in controlled quantum systems such as cold atoms or trapped 
ions has been of great interest. Here we show, both theoretically and computationally, a simplest 
realization of SOC using two-level cold atoms interacting with only one laser beam. The under- 
lying mechanism is based upon the non-adiabatic nature of laser-atom interaction, with the Rabi 
frequency and atom's kinetic energy being comparable to each other. We use the Zitterbewegung 
(ZB) oscillation to further illustrate the effects of the synthesized SOC on the quantum dynamics 
of the two-level cold atoms. We expect our proposal to be of experimental interest in quantum 
simulation of SOC-related physics. 

PACS numbers: 03.75.-b, 32.80.Qk, 71.70.Ej, 37.10. Vz 



I. INTRODUCTION 

Artificial or synthetic SOC in a number of quantum 
systems has been realized. Examples include trapped-ion 
systems [l], 0] .band electrons in graphene Q , cavity elec- 
trodynamics (4, macroscopic sonic crystals and photonic 
superlattices [5(, as well as the more controllable ultra- 
cold atoms in designed laser fields @-|{|. In particular, 
one intensively studied model for implementing SOC in 
cold atoms is the so-called "tripod" scheme, where three 
levels with different magnetic quantum numbers are cou- 
pled with a common excited state by three laser beams 
[a-Q. Therein two degenerate dark states are formed, 
and due to the the center of mass motion of the cold 
atoms, the two dark states are coupled together and as 
a result an effective SOC emerges. However, because the 
three levels on the ground state manifold with different 
magnetic quantum numbers are perturbed by magnetic 
field fluctuations differently, this "tripod" scheme faces 
the challenging issue of dephasing and to our knowledge, 
it is yet to be realized in an actual cold-atom experiment. 
As an exciting progress, cold-atom SOC was realized in 
Ref. pLOj] using a different and somewhat complicated 
scheme, where two of three atomic internal states are 
dressed by a pair of Raman laser beams in the presence 
of a biased magnetic field. 

In this paper we propose the simplest scheme to date 
for the realization of cold-atom SOC, where a two-level 
atom interacts with one laser beam only. In many pre- 
vious proposals for realizing an artificial vector gauge 
potential acting on neutral atoms, the Rabi frequency 
(optical strength) is set to be much larger than other pa- 
rameters of the system (for instance, the atomic kinetic 



energy), such that the optical field is strong enough to 
determine the atomic internal state (dressed state) dur- 
ing the center of mass motion of the cold atom. In other 
words, the internal state of the cold- atom is assumed to 
be adiabatically following the center of mass motion. By 
contrast, here we set the Rabi frequency at a level com- 
parable to cold atom's kinetic energy. This leads to a 
non-adiabatic coupling between the center of mass mo- 
tion and the internal states. Naturally, the non-adiabatic 
coupling depends sensitively on the moving speed of the 
atom and thus creates a type of SOC [ll[ (In Ref. fTTj ) . 
the non-adiabatic transitions or the resulting SOC can 
even induce chaos in the center of mass motion). As 
shown below, the explicit form of the coupling in our 
two-level setup indicates that a simple version of SOC 
can be obtained. 

To elaborate on the dynamical effects of the synthetic 
SOC in our setup, we consider the jittering center of mass 
motion of the cold-atom as an analog of the Zitterbewe- 
gung (ZB) [l2[ of electrons described by Dirac's equa- 
tion. This is of experimental interest because the true 
ZB associated with a real Dirac particle (e.g., electron 
and neutrino) cannot be observed due to its extremely 
large frequency and small amplitude (it does not exist 
even when considering the quantum field theory [l3|). 
Note that artificial ZB has been experimentally imple- 
mented in trapped ions Q , which requires several beams 
of lasers, whereas in our case only one laser beam will 
be needed. This simplification creates more freedom for 
further quantum control purposes. 

This short report is organized as follows. In Section 2 
we shall derive synthetic SOC in a two-level cold atom 
system, of which the two levels are coupled by one plane- 
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wave laser field. In Section 3 we derive the ZB oscilla- 
tion due to SOC. To confirm our theoretical considera- 
tions and to illustrate that the SOC is nontrivial, we also 
show computational results of ZB. Section 4 concludes 
this work. 



II. SIMPLE THEORY 

Consider a two-level atom interacting with a plane 
wave laser beam. The total Hamiltonian reads, 

- _ -ft 2 V 2 / e/2 27LQ co$(kz + wt) \ 

2m + \2hn cos(kz + Ljt) -e/2 ) 

(1) 

where is atom's kinetic energy term, m being the 

atom mass; e is the energy bias between the two inter- 
nal levels; k and oj are the wave vector and frequency 
of the laser field; fio defines the Rabi frequency of this 
laser-atom system. One normally assumes oj = ck, with 
c being the speed of light in free space. However, one 
may also modify this relation by tilting the laser beam. 
Without loss of generality, we assume that the laser beam 
travels in z direction. The total wavefunction can be de- 
noted as = 4>(r)(a,b) T , with 0(r) and (a, b) T respec- 
tively associated with the spatial and internal parts of 
the total wavefunction. The total wavefunction satisfies 
the Schrodinger equation ih-^\ip) = H\ip)- It should be 
noted here that analogous laser-atom Hamiltonian ([T]) 
was studied in the past. For example, it was used to 
study the energy bands of cold atoms in optical lattice 
(e.g., band structure, spontaneous emission of a two- level 
atom, atom's effective mass) (T^ . Our emphasis here is 
placed on the plane- wave states of the cold atom. Physi- 
cally, a wavepacket with sufficient width in the coordinate 
space will be describable by a plane wave. 

We assume that the energy bias e between the two in- 
ternal levels is much larger than the Rabi frequency and 
that the laser frequency is on resonance with the atomic 
transition. We can then exploit the well-known rotating 
wave approximation (RWA) (the validity of this approx- 
imation will be checked in our computational studies be- 
low). That is, under the resonance condition oj = e, 
we make a wavefunction transformation from to 
k//)=<Kr)(a',6') T , 

a' = a cxp(iojt / 2) , b' = 6cxp(— iwt/2). (2) 

In the (a', b') representation, the Hamiltonian (after 
dropping the counter-rotating terms) becomes 

~ -fi 2 V 2 ( hfl Q exp{-ikz) \ 

RWA " 2m + \hQoexp(ikz) J' 

(3) 

We remark that the above representation transformation 
([2|) is independent of the spatial coordinate, so the kinetic 
energy term retains its original form. 

The time-independent RWA Hamiltonian in Eq. (J3j) de- 
pends on both momentum operator — ifiV and coordinate 



z. An important observation is that the Hamiltonian de- 
pends on z through special factors, i.e., exp(±ikz). This 
observation leads to the expectation that the Hamilto- 
nian eigenstates should assume a plane wave form. This 
is indeed the case and can be seen more explicitly af- 
ter making another unitary transformation from \t/j') to 
W) = 0(r)(a",6") T , with' 

a" = a! exp(ikz/2),b" = b' exp(-ikz/2). (4) 

In the (a", b") representation, one can prove from 
the Schrodinger equation that, if the spatial part of 
the wavefunction assumes a plane-wave form <^(r) = 
^Tr^-^eiHP^+P^+P^), the effective Hamiltonian can 
be written in a coordinate-independent form, 

Hsoc = { m ^ Pz ) + 2^ + "8m"' (5) 

which already shows the emergence of SOC (i.e., SOC in 
the form of a z p z , where a z is a Pauli matrix). As seen 
above, the energy splitting in the (a",b") T representa- 
tion depends on momentum p z . It is also seen that the 
momenta in x and y directions do not contribute to SOC 
and that the term is merely a constant. For this 
reason we take p x = p y = for convenience and drop 
the constant. The final one dimensional (ID) effective 
Hamiltonian becomes 

Thus with only one laser beam, ID SOC (SOC asso- 
ciated with one component of momentum) in two-level 
atoms can be realized. Though the well-known "tripod" 
scheme using four-level atoms and three laser beams can 
realize 2D SOC (SOC associated with two components of 
momentum, e.g., Rashba or Drcsselhaus type) instead of 
ID SOC, we note that ID SOC is sufficient for many in- 
teresting purposes and sometimes ID SOC is even more 
advantageous. We also emphasize that in previous ap- 
proaches to synthetic gauge fields, the Rabi frequency is 
always taken to be much larger than atom's kinetic en- 
ergy @ Ell i whereas here the Rabi frequency Slo should 
be comparable to atom's kinetic energy. 

It is also interesting to compare our proposal with the 
realized mechanism in the experiment in Ref. f(Io|). SOC 
in both cases can be interpreted as an outcome of the 
non-commutability between local spin rotation and the 
center of mass motion of the cold atom. The two-level 
set up proposed here is however more straightforward. 
The main assumption made here is that the lifetime of 
the atom's bare excited state is long enough such that 
the spontaneous decay is negligible during a time scale 
of interest to observe the SOC effect. Fortunately, this 
assumption can be fulfilled in real experiments. For ex- 
ample, ytterbium and alkaline-earth atoms have a spin- 
singlet ground state and an extremely long-lived spin 
triplet excited state. Indeed, the lifetime of the spin 
triplet excited state for Yb is almost around 20 s [20l |. 
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III. ZB OSCILLATIONS 

ZB oscillations are now one of the most known effects 
of SOC [1 Hfl El 03 • Interestingly, in the case of ID 
SOC, the damping of ZB can be much slower than that in 
the 2D SOC case @, H, EE E3 ■ Th is difference has been 
explained by analyzing the group velocity of time evolv- 
ing wavepackets undergoing ZB oscillatio ns llal . Follow- 
ing the very same procedure as in Refs. Il6l [l7j . we 
first take an initial wavepacket in (a", b") representation 
as (taking p z = for simplicity), 



llii/ 
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■ 9(z) 7i 



(7) 



with g(z) being a Gaussian center on z = 0. The average 
position of the subsequent matter wave can be analyti- 
cally derived as 



hk 



l Hi 



ze 



-iHi 



sin 2 (f2 t), 



(8) 



where e" 



iHi 



is the time evolution operator. 



In deriving Eq. (|8|), we have made the wide- wavepacket 
approximation as in Ref. 0, IH, HH, [13], i.e., the width 
of Gaussian g(z) in coordinate space Sz is large enough 
such that the width of the corresponding wavepacket in 
momentum space 5p z (~ h/Sz) satisfies 



kSp z 



(9) 



The consequence of the wide-wavepacket approximation 
is twofold. First, according to the Hamiltonian ([6]) the 
spinor will "rotate" around the effective "magnetic field" 
proportional to Hfloi: — j^p z z- The wavepacket in mo- 
mentum space 5p z is very small, so we can take the ap- 
proximation that all the effective "magnetic fields" are 
the same in magnitude (but they are different in direc- 
tions). As a result, the same initial spinor (I,— i) T /\/2 
rotates about different axis with the same frequency 
0, E3 • The evolution for the spatial degree of freedom 
and that for the internal state are thus coupled together. 
Second, according to Eq. ©, the wide-wavepacket ap- 
proximation determines that the amplitude of ZB is much 
smaller than the width of the atomic wavepacket 



hk 



< Sz. 



(10) 



The large ZB amplitude can only arise when Sp z becomes 
larger, or ~ Sl . However, in that case the rotational 
frequencies of spinors are different for different magni- 
tudes of "magnetic field" and Eq. ([8]) does not hold any- 
more, giving rise to non- negligible deforming and fast 
decaying of ZB [7] . 

The small ZB amplitudes are not because of the par- 
ticular realization of ZB in our scheme: it actually stems 



from the quantum nature of atom's center of mass mo- 
tion 0, H, [l6[ (only if the center of mass motion is treated 
classically can the ZB amplitude tuned to be arbitrarily 
large). To enhance the ZB oscillations, one may further 
introduce a driven term in SOC, where the amplitude of 
ZB can be increased or even converted to directed mo- 
tion [l7|. According to Ref. [IT}, we may modulate the 
Rabi frequency periodically to enhance the ZB oscilla- 
tions. We do not pursue such details here. 

We next discuss the (a", b") representation used in our 
theory. This dressed-state representation is related to the 
original bare basis states by two unitary transformations, 
one connected with the optical frequency and the other 
one connected with the laser wave vector. Note that 
both transformations ([2]) and (TJJ do not affect the average 
position of the moving cold atom, so the predicted ZB 
effect occurs in all the three representations. In order to 
make clear the actual physical picture, we now rewrite 
the initial state ([7]) in the original representation when 
t = 0, 



'■'■) 
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1 / exp(-ikz/2) 
y/2 \ -iexp(ikz/2) 



(11) 



which marks the superposition of two wavepackets with 
opposite central wave-vectors. In a lab-frame free space, 
the two wavepackets possess opposite group velocities 
^fhk/m. As a result, without a laser field the initial 
atomic wavepacket will split into two parts. Interestingly, 
in the presence of the laser field, the initial state (fTTj) 
does not split but undergoes ZB oscillation because the 
atomic Hamiltonian is modified by the laser field. The 
initial state preparation here is also simpler than that 
required by early proposals for observing ZB oscillations 
based on the "tripod-scheme [H, [l9| . 

Our theoretical prediction [Eq. (|SJ)] is also verified by 
fully numerical simulations starting from the very initial 
Hamiltonian (JT|) without any approximation. The well- 
known split-operator technique is used in our computa- 
tional studies. One computational result is shown in Fig. 
1, which clearly shows that (i) the ZB period is around 
7r in units of l/f2o (the corresponding frequency is 2 in 
units of fio) and (ii) the ZB amplitude is around 0.5 in 
units of hk/mflo, in full agreement with the theoretical 
result of Eq. ©. 



IV. CONCLUDING REMARKS 

To conclude, we have proposed a straightforward im- 
plementation of SOC using two-level atoms interacting 
with one laser beam. As indicated by our analysis above, 
the SOC coupling strength is easily tunable because it 
is proportional to the effective wave vector of the laser 
beam along the direction of center of mass motion (effec- 
tive wave vector changes when the laser beam is tilted). 
Further, the magnitude of the a x term, which can be 
regarded as the "mass" term of a Dirac particle [after 
neglecting the p\ term in Eq. (6)], can also be adjusted 
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FIG. 1: Numerical result of ZB oscillation with the initial 
state taken as 0. The standard split-operator technique is 
used to simulate the dynamics based on the original Hamilto- 
nian {1} without making any approximation, with e = oj = 50 
(in units of £Iq). In the simulation, the width of the Gaussian 
wavepacket is taken as Sz = 30, which almost keeps a con- 
stant during the shown time scale. Throughout t is in units 
of 1/fio, z is m units of hk/mQo- The numerical results fully 
agree with our theoretical analysis. 



by varying the Rabi frequency, either spatially or tempo- 
rally. 



This work indicates that non-adiabatic coupling in 
laser-atom interaction can be very useful in quantum 
simulation. Indeed, the central idea behind our pro- 
posal is to directly exploit non-adiabatic effects in the 
coupling between the center of mass motion and the in- 
ternal states of long-lived two-level atoms. Our scheme 
is hoped to further simplify experimental studies of SOC 
in the cold-atom context. We have also used the ZB 
oscillations to elaborate on how SOC affects the center 
of mass motion in our setup. It is found that both the 
amplitude and frequency of ZB depend on the Rabi fre- 
quency, which suggests a more flexible control over ZB. 
If we set m ~ 10 _25 Kg, k ~ 10 6 m _1 , we find that a 
ZB amplitude 10 -8 m is obtained if the Rabi frequency is 
taken as ilo ~ 10 5 Hz. A much larger ZB amplitude, say, 
1mm, is also possible if the Rabi frequency is decreased 
to fl ~ 10°Hz, with a ZB frequency about 10°Hz. Both 
the simplicity and flexibility of our two-level-atom scheme 
should be of interest to ongoing studies of spin-orbit cou- 
pling related physics. 
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